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ABSTRACT 


The positions of the ground track of the aircraft and of the aircraft 
sensor footprints, in particular the metric camera and the radar scatter- 
ometer on the C-130 aircraft, are estimated by a program called ACTRK. 

The program uses the altitude, speed, and attitude Information contained 
in the radar scatterometer data files to calculate the positions. The 
report documents the ACTRK program. 
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AND OF THE AIRCRAFT SENSOR FOOTPRINTS 
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Houston, Texas 77058 


1.0 INTRODUCTION 

The purpose of the Aircraft Tracks (ACTRK) program, written In FORTRAN, 

Is to calculate the positions of the C-130 aircraft and of the centers of 
sensor footprints (ground resolution elements or fields of view) during a 
particular data run. The calculations are possible due to frequent (about 
every 0.65 sec) information provided in the data files on the altitude, speed, 
and attitude of the aircraft. In the ACTRK program, it is assumed that the 
aircraft was flying over level terrain, that no errors exist in the input 
data, and that changes in aircraft orientation ^d speed were gradual through- 
out the flight. These assumptions were valid for the agricultural scenes 
addressed by the program. 

2.0 INPUT PARAMETERS 


The input parameters (from the computer compatible tape or disk file) 
are defined as follows: 

1, Roll angle (ROLLM - the angle between the line passing through the 
aircraft wing tips and the horizontal plane. The roll angle is positive when 
the left wing is above the horizontal plane. The units are degrees. 


^The name of the parameter in the FORTRAN program is given in all capitals. 
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2. Pitch angle (PITCH) - the angle between the line passing through the 
fuselage of the aircraft and the horizontal plane. The pitch angle is positive 
when the nose of the aircraft is above the horizontal plane. The units are 
degrees. 

3. Heading angle (AHEAD) - the azimuthal orientation of the line passing 
through the fuselage of the aircraft as f ‘en from above relative to true north. 
The angle is measured clockwise from true lorth. The units are degrees. 

4. Drift angle (DRIFT) - the differei e between the heading angle and 
the direction of the actual ground track. The drift angle is positive when 
the ground track is clockwise from the heading angle. The units are degrees. 

5. Ground speed (GS) - the speed of the aircraft with respect to the 
ground. Ground speed is calculated by the onboard navigation system computer 
from the air speed and the wind speed vectors. The units are knots. 

6. Radar altitude (ALT) - the vertical distance between the aircraft 
and the ground. The units are feet. 

The above parameters are given on the first of the four data files versus 
elapsed time (ET) from the time of the first record. Also, the actual times 
(Greenwich Meridian Time, GMT) of the acquisition of these data are given on 
each record in hours (NHR), minutes (NMIN), and seconds (SEC). These time 
values are used to compute the time in seconds (TIME). The units for time 
(ET and TIME) are seconds. 

3.0 LOGIC OF THE PROGRAM 

If one follows the logic of the ACTRK program, it is divided into several 
parts as follows. 
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Dlmenslons of arrays are defined at Che beginning. Allowances are made 
to handle up to 300 records (data records) for each aircraft run. This can 
be increased if necessary. 

The program requests an irput from the terminal for the orientation of 
the flight line (FLIGHT) that was being flown (the intended flight line 
direction measured clockwise from true north as viewed from above). Be sure 
Co Include a decimal point in Che input from the terminal since the read 
format is floating point. The units are degrees. 

The flrsL three records of the file are read into dummy variables. 

These records are headers for the remaining data records. 

Data records are then read in a loop until the end of file is encountered. 
Within the loop, the time data are converted to seconds (TIME), and the de- 
rived parameters, GTRACK and G/JIMA, are calculated. GTRACK is the ground 
track angle (in degrees from true north in a manner similar to that of AHEAD). 
GTRACK is given by 

GTRACK - AHEAD + DRIFT (1) 

GAMMA is the angle between the flight line direction (FLIGHT) and the ground 
track (GTRACK) and is given by 

GAMMA - FLIGHT - GTRACK (2) 

The ground speed (GS) is converted from knots to feet per second and re- 
stored in tne same variable name. 

After the Inputting of Che data records is completed, Che number of data 
records, ISTOP, is calculated. Then, two more inputs are requested from the 
terminal. One is the absolute position (XSTAR7) of the aircraft at the time 
of Che first data record as measured down the flight line (X-axis). The 
other is the absolute position (YSTART) of the aircraft at the time of the 





first data record as measured to the left of the flight line (Y-axis) as one 
faces the flight line direction. Be sure to Include decimal points in these 
parameters entered from the terminal since the format is floating point. The 
Incl kSlon of XSTARl and YST/.PT allows one to locate the subsequent positions 
in an absolute frame of reference. The units are feet. 

The program proceeds by going through a loop in which the positlcns 
(X,Y) of the aircraft are calculated. A simple dead reckoning procedure is 
used. The aircraft is vectored along at the average ground speed (ari‘<’.hmetic 
average of two adjacent data record values) in the average direction anjle, 
GAMMA, with respect to the intended flight line, FLIGHT. The time of travel 
is the difference in the elapsed times (ET) between data records. Elapsed 
time is used rather than absolute time (TIME) since ET is given to an accuracy 
of hundredths of a second and TIME is given only to an accuracy of tenths of 
a second. The results of the aircraft position calculations are stored in 
the X and Y arrays. The units are feet. The equations for X and Y are as 
follows : 

^i+1 “ Xi + (ETi+i - ETi) * Ci * COS( GAMMA) (3) 

''^i+l " ^i (ETi+i - ETj) * GS * SIN(GAMMA) (4) 

Given the aircraft positions, X and Y, rhc program proceeds with the 

calculation of the positions of the photo centers and the centers or the 
radar scatterometer footprints (i.e., the areas from which returned and 
filtered radar scatterometer signals come). Another angle, EPI, is calcu- 
lated in this procedure. EPI is the difference between the 'ntended flight 
line direction (FLIGHT) and the heading angle (AHEAD) of the aircraft. EPI 
is needed due to the fact that the f^n-shaped antenna of the radar scatter- 
ometer is aligned with the axis of the aircraft and is directed toward the 
aft of the aircraft. Aircraft roll affects the location of the intersection 
of the antenna pattern and the ground. It is assumed chat Che metric camera 
is hard mounted to the aircraft (i.e., it rolls with the aircraft as the air- 
craft rolls). In the case of the radar scatterometer, changes ii. pitch angle 
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(PITCH) do not aff«ct the position of the radar scatterometer footprint since 
the Incidence angle has been specified according to the Doppler freouency 
shift filter used in processing the data. Changes in pitch angle do affect 
the position of the cauera phot^ center. Changes In the roll angle do affect 
both the camera and radar acatterometet footprint positions. Of course, 
changes In altitude affect these posltlonit also. In the last computational 
loop of the program, the equations used to predict the effects of changes in 
pitch, roll, and drift on these positions are applied. See the attached 
figures (Figures 2 and 3) for the equations. 

Ihe last loop simply writes the results of the calculations out to a 
file for subseuqent use. 

The last plot given In this note shows the results of the application of 
ACTRK to an actual aircraft flight run over Webstei County, Iowa, on August 19, 
1980. Note the different tracks of the photo center, aircraft ground track, 
and radar scatterometer footprints. In some Instances, the lateral displace- 
nh>nt between pairs of tracks Is several hundred feet. This could result In 
some footprints lying outside a particular field of Interest at some or all 
angles of viewing. It should be noted. In particular, that the use of photog^ 
raphy to locate the aircraft and the positions of radar scatterometer foot- 
prints could result in some data being included wrongly in a data set that is 
supposed to be taken over a given field. 

4.0 CONCLUSION 


The proc'sm, ACTRX, is explained in this short document. It produces 
useful estimates of the positions of the aircraft, photo centers, and radar 
scatterometer footprints where level terrain has been overflown and where 
changes in the altitude, speed, and attitude of the aircraft are slow enough 
to be well tracked by the data given In the first data file of the set of 
four data files produced by the NASA Johnson Space Center before 1982. 
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If significantly unlevel terrain should be encountered, the user should 
not use the ACTRK program. The distance from the aircraft position to the sen- 
sor footprint position would be in error. Also, the range from the aircraft 
to the sensor footprint would be in error. Finally, the local incidence angle 
would not be equal to the sensor look angle. It should be possible to alter 
the program such that the elevation of the surface along the flight line can 
be estimated on one pass and used in antoher pass to make these estimates 
more accurately than dues the existing ACTRK program. 

To use the calculated aircraft and aircraft censor positions with the 
measured data, the author makes the following comments. One can use the air- 
craft positions versus time of flight (TIME) or elapsel time (ET) as given as 
a basis for other comparisons. The sensors, however, acquire data over an 
area that is not, in general, directly below the aircraft. For example, a*: 
nominal ground speeds (160 knots), at an altitude of 460 meters (1500 feet), 
and with a pitch angle of one degree, the center of the photos leads the air- 
craft position by 8 meters (26 feet) or 0.1 seconds. The position of the 
radar scatterometer footprint taken at a sensor look anglo of 15 degrees lags 
the aircraft position by 123 m (402 fc.) or 1.5 sec. At a sensor look angle 
of 50 degrees, the lag is 548 m (1788 ft) or 6.6 sec. Reportedly, in the 
case of the radar scatterometer data, the measurements have bee.i adjusted 
such that the radar backscattering c^'afflcient value reported at a given time 
is of the area under the aircraft vor to the left of right of the aircraft 
when nonzero roll and drift angles existed) at the time reported. This is 
the so-called nadir time correction. Nevertheless, the data are actually 
acquired several seconds later (1.5 and 6.6 seconds, respectively, for 15 
and 50 degree data) when the aircraft wo<? in a different attitude than at the 
time of the overflight of the area in question. The author has noted, also, 
that the apparent time that the airplane crosses a field boundary according 
to the photography can be as much as 1.5 seconds before that boundary is 
crossed by the radar scatterometer even though the data have been nadir time 
corrected. The explanation of this observation has not been found by the 
author. 
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FIGURE 1. GEOMETRY OF AIRCRAFT ATTITUDE 

T - GAMMA ANGLE 
c ■ EPI ANGLE 
6 - DRIFT ANGLE 
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AIRCRAFT 

POSITION 


FIGURE 2. RELATIONSHIP BETWEEN PHOTO CENTER AND AIRCRAFT 
POSITION 

A - ALTITUDE (FEET) 
p - PITCH ANGLE (DEG) 
r - ROLL ANGLE (DEG) 
a ■ A tan p 
b « A tan r 

^p * ^N e - b cos E 

Yp * Yfj *• a sin e + b cos e 
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FIGURE 3. RELATIONSHIP BETWEEN RADAR FOOTPRINT POSITION AND 
AIRCRAFT POSITION 

A - ALTITUDE (FEET) 
r - ROLL ANGLE (DEG) 

G - RADAR LOOK ANGLE (DEG) 
b ■ A tan r 
c » A tan 0 

d *y 

Xr * % ~ cos e - b sin e 
Yr * Yn ” d sin e + b cos e 
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T“tE OBJECTIVE OF TdlS P^^uOPav is 
AlKCPAFT (aPPABENI Oh?0UN() TPaCK) 
OF THE wtTPlC camera, and of Th£ 
SPECIFIC angles. TfiE INPUT DATA 
speed. ANn altitude INFuKMAriON 


ACT 

TO calculate The POSITIONS OF ThEaCT 
. OF TmE principal POiM (CEMtP) ACT 
PAOAK SCATTEPO**ETEr* FOOTPRINTS aTaCT 

consist o^ aipcraft attitude. act 

LOCAFEu IN_TPE first FILE OF ACT 


adar scattepometep data ppqvideo'ow computer compatible Tapes bt 

ThE NASA JOHNSON space CEnTE^ oEPDkE iSPE. 


THE PROGRAM' WAS WRITTEN BY DP. 
TEaAS. on StPTEMotP 7. 19H3. 


JACK F. Paris. NASA/JSC. HOUSTON. 
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UOOlO 
OOOiiO 
OOUSO 
UOOuU 
oOOSO 
OOOoO 
GOO 70 
OOOBO 
UOOVO 

uoloo 

oOlIO 

00120 

00130 
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DEFINITION OF parameters 


parameter 

ET 

TIME 

poll 

PITCH 

DRIFT 

ALT 

G5 

ahead 


flight 

GTRACK 


GAMMA 

tPl 

XSTA«T 


ystart 

X 

Y 

XR. YP 
X15. Y15 

X50. Y50 


definition 

ELAPSED time (SEC) SINCE PlJN START 

CLOCK TIME (SEC) FROM GMT TIME ON DATA RECORDS 

roll angle (DEd) (POSITIVE IS RIGHT wInG DU<«N) 

PiTcn angle (oeg) (Positive is nose up) 

JrlKT ANGLE (OEG) . GROUND TRACK HEAdInG MINUS 

aircraft heading 

altitude of aIpCPAFT (FEET) MEASURED OY RADAR 
altimeter (ALONG VERTICAL TO GROUND) 
jROUND SPEED (KNOTS) CALCDLATEd BY ON-BOARO 
NAVIGATION SYSTEM 

aircraft HEADING (DEo) MEASURED FROM TRUE NOpTh 
(0 - 360 DEG). IS NOT ThE SAME, IN GENERAL. 

AS The ground track 

ThE intended GROUND TRACK OF THE AIRCRAFT OR THE 
DIRECTION (DEG) OF ThE FlIGhT LINE MEaSURED F RQm 
fPUE (NORTH (0 - 360 DEG) 

AIRCRAFT HEADING PLUS DRIFT (OEG). The ACTUAL 
(■‘IRECTION THE AIRCRAFT MAKES RELATIVE TO The 
o?OUNU measured from TRUE NORTH (0 - GbO DEG) 
FLlonT .MINUS GROUND TRACK (dEG) 

FlIoHT mInjS aircraft hEauIivG (DEG) 

ABSOLUTE position OF' ThE AIRCRAFT alOnG ThE 
flight line at the timc of The first uATA 

pECOkD for The RuN (F'fcET) 

ASSULUTE POSITION OF THE AIRCRAFT 
lu The FLI^’Hf line aT ThE TIKE OF 
mECOpO fop ThE PUN (FEeT) 
position OF The aIrCRAFT MEASURED 
fLIomT LInF. (FtET) 

POSITION OF The aircraft measured perpendicular 
to ThE flight line (FEET) 

POSIT UiN OF The center OF THE METRIC PHOTO 
hOSITICN of the 15 OEG RADAR SCaTTEhUkETER 
FOOTPRINT (FEET) 

PoSITlCtN OF T h£ so deg pADar SCATTLpOpETER 
footprint (FEET) 


perpendicular 
The first data 

ALONG The 


(FEET 




DIMENSION ET (30n 

1. drift (300) .GTkA 

2. »1S(300) .YiS(30 


) .mOllO^iO) .pITCh(300).AlT(j 00).GS(300) .AHEAD ( jOO 
Ck (J.'O) .GAMvA ( 300 ) .X (300) . Y ( 300) .XR (30u) . YP ( jou) 
0) ,xou(30u) .Ybu(300) .TIME (JOO) 


Input of flight lINe uI-fCTION from The TEf'MiNAL 


• p 1 T E ( 1 6 « .? 0 
200 format (• FN 
READ( 15. lOy 
loo format (F li--. 


TE-y iNTENUtO aircraft FlIghT LInE HeAOInG (DeG):') 
) Fl K’hT 
2 ) 


convert A':'jL&n From DEj^EES TO pAOIA.o 


Pls 3 . 1 A 15 P 2 
The 15W = 1 5 . 
TmE 50 P = 50 .•» 


Pl/lH 
PI/ Ir 


dummy peao FiPST r-t-'t'L -^eCOPEj (heauERS) OF THE DaTa FIlE 


ACT00160 
ACT00170 
aCTOOIbO 
ACTOOl^O 
ACTU0200 
aCTOOBIO 
ACT00220 
ACTU0230 
aCTOOEmO 
aCTu0250 
ACT00260 
ACTU0270 
ACT002%0 
ACT00290 
aCTOOJOO 
ACT00310 
ACT00320 
ACT00330 
ACTOOJ^O 
ACT00350 
ACT00360 
ACT0037& 
ACTU03S0 
aCT003S»(; 
aCTOOmOO 
ACTOOMIO 
ACTU0‘»2C 
ACT00M3C 
ACT OOMMt 
ACTOOA5C 
aCT 00'»6C 
ACT0047C 
ACTOOAaC 
ACT00A9C 

actoosoo 

aCTOOSIO 
ACTU0520 
ACT00530 
) ACT00540 
aCTOOSSO 
ACT00560 
ACT00570 
ACTOOSdO 
aCTOOSPO 
•ACTOO&OO 
ACT00610 
) ACT00620 
act U0b30 
ACTOObRO 
ACTUObSO 
aCTo O bbO 
ACTu0b70 
ACTOOboO 
ACTuObRO 
ACTuO 700 
ACT v07l 0 
ACT 00720 
ACT u0730 
aCT u07>*0 
aCTw0750 
aCT II 0760 
ACTu0770 

aCTo 07 p 0 
aCTOO /so 
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f-'lLt: ACTf^K FOf^T^AN A EODL / JOHNSON SPACE CEnTES 


ALSO. DJMMY r^cAD F I^ST DATA PfcCOPQ aINCt ET DIFFE/^EnCE d£T«EtN 
Flt^ST and second data kECORO IS AlSO TOO LA.^0E 

PEAO(20.101 ) AA.Bd.CC.uD 
lul format (A5/A5/A5/A3) 

BEGIN LOOP TO READ OaTA PECORUS FROM ThE DISK FILE 
DO 20 1=1. JOO 

»EAD(20. lu 2. E. >10 = 300 )ET ( I ) .NhR.NmIN, SEC.ROLu ( I) .PITO-K I ) .DRIFT ( I ) 
1 .alt ( I ) ,GS ( I ) . AM cAO ( I ) 

102 FORMAKr 7.2.U.I2.1X.I2.1X.F5.2.F12.2.5F10.2) 

CONVERT MP-mIn-SEC TO SECONDS (£.mT) 

TIME ( I ) =3SO.).»FLOAT (NHW) ♦ esO.^-'FLOAT (NMIN) ♦ SEC 
•/RITE INPlT Data TO lINE PRINTER TO VEkIFY C0«RLCT INPUT 

NOTE **•* <j5Er '^ay Eliminate tnIs ste® 

•(RITE (6.5f)0)ET (1) .NMr.nmjn, sEC.kOlL ( I ) »PI Tph ( 1 ) .Or IFT ( 1 ) . alT ( I ) 

1 »GS ( I ) . AHE AL { I ) 

500 format (FB. 2.2 13. Fb. 1 .bFt-.2) 

calculate ground track and gamma angle 

GTRACk ( n =ArlEA0( I ) ♦ urIFT(I) 

6AMMA(I)=FLIGHT - OTRACK(I) 

CONVERT GROUND SPEED TO FEET PER SECOND FROM KNOTS 

US(I) =GS(I)<*60bO./3SOO. 

20 CONTINUE 

ISTOP IS number of data RECORDS (ASSUMED TU BE LESS THAN 300) 

300 IST0P=I- 1 

Enter absolute position of aircraft at timl of first data record 

WRlTEdB.POU 

201 FORmatc enter XSTaRT (FEET):') 

REA0(15.in0)XSTART 

WRITE (lb. 202) 

202 FORMATC EnTFR YSTART (FEET):') 

READ (15.100) YSTaRT 

X (1)=XSTAPT 
Y ( 1) =YST APT 

calculate Increment of time, average ground speed, average gamma 
angle, and nE-T aircraft POSITIuN 

00 21 1=2. ISTOP 

lMi=I - 1 

0ELTAT=£T(I) - ET(IMI) 

GS A VG = 0 . ( GS ( I ) ♦ OS(lMl)) 

GAAVU=0 , 5* (gamma ( I ) ♦ GAMMAdMl)) 

GAAVf 5 P = GAAVG^Pl/lo'i. 

DELX = GSft VG^UElTAT^'COS (GAAVGR) 

DELY=GSAVr.<»uELTAT<»5lN(GAflVGR) 

X ( I) =X ( IMl ) * DELA 

Y(I)=Y(IM1) . DELY 

21 CONTINul 

CALCULATE PROTO POSITION AnO rAuAP bCATTERU^'ETER FOOTPRINT 
POSITIONS (15 AND SO UEGREES) 

00 30 I=l.l>TOP 
EPI=FlIGHT - ArtEAUd) 
tPlR=ERI'^PI/l'<0. 

SEPIR=SIN(£rIR) 

CEPlP=C0S (EPIP) 

P®=PITCn ( d RPI/IRO. 

RPsROLL ( I )MP i/lB0. 

A=ALT ( I ) *TAN (Dr) 

B = ALT (I)»TAf>.(RR) 

C15=ALT ( I ) MTAN ( TRt 1 3- ) 

CSQ=ALT(I)«Tan( TrESOR ) 

D15 = SC'RT (Clb^ClS - 
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ACT01300 
ACT01310 
ACT01320 
ACT01330 
ACT01340 
aCTo 1 35o 
ACT01360 
ACTO1370 
ACTol3bO 
ACT0I3R0 
ACT 01400 
ACI 01410 
ACT 0 I 420 
ACIyl 430 
A C T 0 i *.4 0 
aCTo 14'50 
AC I Ol4bO 
ACT01470 
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AC T 0 1 <.9 0 

ACT0I500 
act 0 IbiO 
aC T u 1 b20 
act 01530 
aCTO I D40 
ACT01550 
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ACTuib70 
aC TO 1 5bu 



nno 


ORIGINAL PAGE tS 
OF POOR QUALITY 


fill*. ACTRK F0»T*^AN a LuOl / JOhNSOa) space center 


DSOsSORT (C5O*C50 - ^*>3) 
XD ( 1 ) =X ( I ) ♦ a*CeP IR - 

YP ( 1 ) =Y ( 1 ) ♦ A* 5 £PIR ♦ 

X15(I)SX<I) - J15»CEPIR 
X50 ( 1 ) =X ( I ) - J50*CEJ1R 
Y15(I)=Y (I) - DIS^SEPIR 
YSO ( I ) =Y ( I ) - DS0*»S^PIR 
3() continue 


- 5 *SEPIh 

- B»Sc.PIR 

- P* 5 EpIR 

♦ 0 RCEPH 

♦ H«^CEPIR 


*<*^1TE results of calculations TO DlS^ E IlE 
00 22 I=l.lSTop 

'^-'ITE { 3u» 203) X 1 1 ) , Y ( I ) ,aP n ) * >P ( I ) »X15 ( I ) . Y 15 ( I ) , I I ) , Y^>0( 1) 
2U3 format ( U.hF 9.0) 

22 continue 

STOP 

End 


ACToISrO 

actlI&oo 

ACTui'slU 
ACT01620 
ACT0163U 
aCToIo**© 
ACTOlGSO 
ACTO 16b0 
ACTUlfeTO 
ACTOiobO 
ACT Olb>yO 
ACTOl fOO 
ACT01710 
ACT01720 
ACTOl 73C 
ACTOl 7**u 
ACTOl /SO 



